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Q,' ABSTRACT 

^ : 

^ , We reply to the recent criticism by Garriga and Tanaka of our proposal that 

quantum gravitational loop corrections may lead to a secular screening of the 
^ I effective cosmological constant. Their argument rests upon a renormalization 

t:;-]- ■ scheme in which the composite operator {Ry/^ — 4Ay^)i.en is defined to 

O ■ be the trace of the renormalized field equations. Although this is a peculiar 

prescription, we show that it does not preclude secular screening. Moreover, 
we show that a constant Ricci scalar does not even classically imply a con- 
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O I stant expansion rate. Other important points are: (1) the quantity i?ren of 

Garriga and Tanaka is neither a properly defined composite operator, nor is it 
constant; (2) gauge dependence does not render a Green's function devoid of 
. physical content; (3) scalar models on a non-dynamical de Sitter background 

rS I (for which there is no gauge issue) can induce arbitrarily large secular con- 

I tributions to the stress tensor; (4) the same secular corrections appear in 

observable quantities in quantum gravity; and (5) the prospects seem good 
for deriving a simple stochastic formulation of quantum gravity in which the 
leading secular effects can be summed and for which the expectation values 
of even complicated, gauge invariant operators can be computed at leading 
order. 
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1 Introduction 



Some years ago we proposed that the continual production of infrared gravi- 
tons induces a sort of quantum friction which gradually slows inflation [1]. 
When inflation rips a pair of infrared gravitons from the vacuum, which 
is a 1-loop effect, they induce a gravitational potential that contributes to 
the energy density at next order. This potential remains imprinted on the 
spacetime even after the pair has been pulled out of causal contact. The 
number density of pairs remains constant because the vast expansion of the 
spatial volume cancels the continual creation of pairs. However, the induced 
gravitational potential grows without bound. Taking account of the num- 
ber of infrared gravitons inside the past light-cone of a local observer gives 
$ ~ —GH^ ■ Ht [2j. The effect of this potential must be to slow inflation 
because gravity is attractive. Hence the induced energy density of interaction 
is negative, p ~ —GH^ ■ Ht. Because gravity is weak, this energy density 
grows very slowly, in other words, |p| ^ H\p\ for Ht ^ 1. Combining that 
fact with stress-energy conservation, p = —3H{p + p), implies the induced 
pressure is p ~ —p. Hence one gets a growing, negative vacuum energy. 

It is tempting to speculate that this mechanism might simultaneously ex- 
plain why the observed cosmological constant is so much smaller than the 
natural scales of particle physics, and also provide a model of inflation which 
has no fundamental scalar. The idea is that primordial inflation would start 
because the bare cosmological constant is positive and of GUT scale. This 
avoids the problem of needing a scalar inflaton to be unnaturally homoge- 
neous over a super-Hubble volume Inflation would persist for a long time 
because gravity is a weak interaction, thereby dispensing with the need for a 
shallow potential. And inflation would be brought to an end by the gradual 
accumulation of negative vacuum energy. 

Explicit computations in this scheme are challenging because the first 
effect should be at two loops and because the putative effect would be non- 
perturbatively strong during the current epoch. Ford was early able to show 
that there is no secular deviation from de Sitter background at one loop [3], 
in agreement with the mechanism. This 1-loop result has been confirmed 
by two other groups [5l[6]. A year-long computation of the graviton 1-point 
function revealed secular slowing at the expected 2- loop order [7j|3 Because 
back-reaction should be small until the last few e-foldings, it was possible 

^This work has not been confirmed, nor re-done using dimensional regularization [6]. 
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to use the 2-loop result to estimate the scalar and tensor power spectra [S]. 
Some thought has also been given to how the mechanism might operate after 
it becomes non-perturbatively strong at the end of inflation [9]. 

Our proposal has recently been criticized by Garriga and Tanaka flU\ . 
They raise two objections: 

1. The graviton 1-point function cannot show slowing, or anything else, 
because it is gauge dependent; and 

2. There can be no slowing of inflation in pure quantum gravity because 
the field equations set the renormalized Ricci scalar to a constant. 

In section 2 we point out that physical information can reside in gauge depen- 
dent quantities, that the secular growth we found was not injected through 
our gauge fixing functional, that the only invariant quantities so far checked 
in quantum gravity show the same type of secular dependence, that this sec- 
ular dependence also affects the vacuum energy of scalar models for which 
there is no gauge issue, and that it could only drop out of the quantum grav- 
itational vacuum energy through an infinite series of cancellations. Section 
3 contains two important points: 

1. Garriga and Tanaka did not show that the renormalized Ricci scalar is 
constant and the relation they actually derived is completely consistent 
with secular slowing; and 

2. They did not define the Ricci scalar as a properly renormalized com- 
posite operator but rather the Ricci scalar times the volume element. 

We also comment on some peculiarities of their renormalization scheme. Sec- 
tion 4 makes the short but crucial point that a constant Ricci scalar does 
not even classically imply expansion. In Section 5 we review recent progress 
in developing a computational formalism which is powerful enough to evolve 
into the non-perturbative regime. 

2 Gauge Dependence 

The fact that the expectation value of the metric cannot be formed into a 
gauge invariant quantity in the same way as a classical metric was originally 
pointed out by Unruh [llj. We agree with him, and two years ago we ad- 
vocated that the quantum gravitational back-reaction on A-driven inflation 



2 



should be quantified with an invariant operator which measures the geodesic 
deviation [12] . The expectation value of this operator has not yet been com- 
puted at 2-loop order owing to its complexity. It is worth noting that the 
difficulty of working with gauge invariant observables has led Losic and Un- 
ruh as well to the consideration of gauge dependent quantities in their own 
interesting study of back-reaction [T3] . 

• Gauge Dependence Does Not Automatically Imply Lack Of Physical 
Content: A good example is provided by the IPI functions of the Standard 
Model. These IPI functions are certainly gauge dependent, yet they can be 
assembled to give the gauge independent S-matrix, which has of course been 
subject to impressive experimental verification. So the gauge-dependent IPI 
functions of the Standard Model must contain valid, physical and gauge- 
independent information along with some unphysical and gauge-dependent 
behaviour. 

An especially noteworthy example is the gauge dependence of effective 
potentials [H]. In the days when people were studying the possibility for 
loop corrections to stabilize Kaluza-Klein compactifications P^j this gauge 
dependence was regarded with great suspicion and even despair. The ques- 
tion was asked, "which result should be trusted if the effective potential 
shows a non-trivial minimum in one gauge but not in another?" It was even 
believed that a "field-space metric" could be identified that would allow con- 
struction of a gauge independent (and also field re-definition independent) 
effective action [T7] . 

It was eventually realized that gauge dependence is actually a blessing for 
this problem. What one really wants is to find a stationary point of the full 
effective action; it is only because the effective action is too complicated to 
evaluate for an arbitrary background field that one is reduced to studying the 
effective potential. Because different gauges are related by field-dependent 
gauge transformations, the same finite-parameter family of field configura- 
tions in one gauge can, in another gauge, probe a different direction in the 
full space of fields. Hence the correct conclusion is not that both effective 
potentials are rubbish on account of their gauge dependence, or that one is 
right and the other wrong, but rather that each one represents valid physical 
information about the theory. For a putative solution to be correct it must 
make each of the effective potentials stationary — and even then it might still 
be spurious if there is a non- stationary direction that neither of the restricted 
probes of field space chances to access. 
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The expectation value of the graviton 1-point function in a fixed gauge is 
a perfectly valid measure of quantum corrections to the background in that 
gauge. It might be that the secular screening we found is an artifact of the 
gauge (or even an error in what was a very long and difficult computation) 
but several possibilities can be ruled out. First, the secular dependence we 
found was not introduced through the gauge fixing functional. We found "in- 
frared logarithms" ; these are logarithms of the ratio of the scale factor at the 
observation time to its value at the beginning of infiation. These logarithms 
break the part of the 10-parameter de Sitter group known as dilatations; the 
latter generate the transformation — > k ■ in conformal coordinates. 
However, our gauge condition preserves this symmetry. 

• Infrared Logarithms Do Not Always Drop Out Of Physical And Gauge 
Invariant Quantities: The only invariant so far checked is the infiationary 
power spectrum and Weinberg has shown that infrared logarithms do correct 
this |18[[T9]. Because the logarithms in that case are of the ratio of the current 
scale factor to the scale factor at horizon crossing, their enhancement is at 
most about ln{a/a^c) ~ 100 for any perturbation whose spatial variation 
we can resolve. This is not enough to overcome the small loop counting 
parameter of GH^ < 10~^^. However, there can be arbitrarily large infrared 
logarithms correcting things we perceive as spatially constant such as the 
vacuum energy. 

It might still be that the vacuum energy is somehow protected from 
the secular effects which contaminate IPI functions |71 [20], the quantum- 
corrected mode functions [21] and the power spectrum [181 [19]. But that is 
certainly not true for scalar models on non-dynamical, de Sitter background. 
In that case the calculations are vastly simpler, and there is no gauge issue 
to frustrate drawing obvious physical conclusions. It suffices to compute the 
expectation value of the stress tensor. For example, a fully renormalized, 2- 
loop computation of the stress tensor of a massless, minimally coupled scalar 
with a ip'^ self-interaction which is released in Bunch-Davies vacuum at t = 
(with scale factor a{t) = e^*) reveals the following secular growth for the 
induced energy density and pressure [22]: 




(1) 



Pit) 



-{-lln'{a)-^Ha)} + 0{X'). 



(2) 
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At leading logarithm order this result is in perfect agreement with the non- 
perturbative analysis of Yokoyama and Starobinskii It also has the 

transparent physical interpretation that the inflationary production of scalars 
increases the scalar field strength, which drives the scalar up its potential 
and thereby induces a growing vacuum energy. 

One may dismiss the ip^ result on the grounds that it shows an increase 
in the vacuum energy, whereas we claim a decrease for quantum gravity. But 
it is simple to find scalar models whose induced vacuum energy is negative. 
Consider scalar quantum electrodynamics, for which the analogous 2-loop 
results are I25l ES : 



pit) = - — — X 



-^ln(a) 
^ In (a) 



+ O(e^) , (3) 
+ 0{e'). (4) 



The physical interpretation seems to be that the inflationary production of 
charged scalars polarizes the vacuum [271 [28] , which lowers the scalar energy. 

Again one may dismiss the scalar QED results on the grounds that 
they represent a transient effect that approaches a constant after a non- 
perturbatively large time At ~ l/{e^H) [26]. The specious nature of this 
argument can be seen from a massless fermion which is Yukawa coupled to a 
massless, minimally coupled scalar on non-dynamical de Sitter background. 
For that model the induced vacuum energy falls without bound [29]. This 
also has a transparent physical interpretation: the inflationary production 
of scalars increases the scalar field strength, which induces a fermion mass, 
whose effect is to decrease the vacuum energy. Had gravity been dynam- 
ical in this model, the resulting secular back-reaction would end inflation. 
That the universe subsequently decays to a Big Rip singularity does not al- 
ter the model's demonstration that secular back-reaction can accumulate to 
dominate late time cosmology. 

Finally, it should be noted that the absence of secular quantum gravita- 
tional back-reaction would not only require that a single infrared logarithm 
drops out at two loops. The number of infrared logarithms grows with each 
loop in a way that can be predicted from the number of undifferentiated 
massless, minimally coupled scalars or graviton fields in the basic interac- 
tion vertex. The number of extra infrared logarithms for each extra coupling 
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constant in Xip'^, Yukawa theory and scalar QED are [2^] : 

These expectations are verified in a large number of explicit 1-loop, 2-loop, 
and even 3-loop computations [30]. They follow as well from the stochastic 
analysis of Starobinskii and Yokoyama [231 l2ll [26] . 

The basic interaction of quantum gravity takes the form K'^h"'dhdha^~'^ 
|32j . where k"^ = IGttG is the loop counting parameter. The number of extra 
infrared logarithms for each extra factor of GH^ is |26] : 

K^'h^dhdha^-'^ =^ In (a) for every GH"^ , (8) 

So the absence of secular back-reaction in quantum gravity would require that 
the escalating series of infrared logarithms which do appear in IPI functions 
and physical quantities somehow contrives to cancel out of the expansion 
rate. We do not regard it as reasonable to suppose that this happens. Of 
course Garriga and Tanaka claim to have proven that it does happen, and to 
all orders. We turn now to an explanation of why their proof is not correct. 

3 Renormalizing the Ricci Scalar 

The basic argument of Garriga and Tanaka consists of a peculiar scheme for 
renormalizing the Ricci scalar so that the following relation applies: 

(i?rcnV^ ) = 4A ((l + 5™l)v^ ) • (9) 

We shall presently describe some technical problems with this scheme but 
the most significant point of this section is that Q does not imply a constant 
Ricci scalar. We shall demonstrate that Q is completely consistent with the 
relaxation mechanism of Section 1. Moreover, R^^n as defined by Garriga and 
Tanaka is not a properly renormalized composite operator because it fails to 
produce finite results when inserted in IPI functions. 



^Certain IPI functions have counterterms that can be used to absorb mfrared loga- 
rithms at low loop orders |31j but factors of ln(a) always show up at higher orders. 



-g =^ In (a) for every 
-g =^ ln(a) for every 
-g ln(a) for every 



A, (5) 
f , (6) 
. (7) 
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• Preliminaries: Many people have concluded that screening is impossible 
in pure quantum gravity because the trace of the classical field equations 
imphes i? = 4A: i 

'^^li^i"-'")-^ ^ 16.G9-^ = (-fl+4A)v/^ = 0. (10) 

The conclusion is premature because quantizing general relativity requires 
adding to the Lagrangian an infinite sequence of ever-higher dimensional, 
BPHZ counterterms: 

AC = a^R^^ + a2C""^^"C,.M-v^ + • • • (11) 



The divergent parts of the Oj's are fixed by the need to cancel primitive 
divergences, but no physical principle seems to fix their finite parts, and 
what they are affects physical predictions. That is one way of expressing the 
problem of quantum gravity. 

Of particular interest to us is the fact that even the lowest order coun- 
terterms alter the relation between R and A implied by the trace of the 
renormalized field equations: @ 

= IGnGg^"^^^^^^^ = (-R + 4A + IGttG x GaiDR + . . )^ . (12) 



One would typically break the a^'s up into divergent and finite parts, and 
render this equation in terms of suitably renormalized composite operators: 

0= (-R^/^+AA^/^) +167rGx6ai finite fn^v^) +0(G^) . (13) 

V / ren ' V / ren 

The Eddington counterterm — the one proportional to ai — was exploited 
by Starobinskii in constructing an early model of what would later be called 
infiation ^34j. It should therefore be obvious that the equations of quantum 
gravity cannot imply the Ricci scalar is constant unless the finite parts of a 
countably infinite number of a^'s are set to zero. 

• The Garriga and Tanaka Argument: They propose to circumvent the 
problem of counterterms by absorbing them into the renormalized Ricci scalar 



•^We first heard the argument from L. Susskind in the mid 1980's [53]. 

^Computing the trace in dimensional regularization would result as well in finite con- 
tributions from the ai and a2 counterterms whose inclusion would only strengthen the 
argument we shall make. 
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i?ren and the renormahzed volume element (1 +S^oi)\/—g- Their equation 
(36) defines the operator whose vanishing in expectation values is claimed to 
preclude secular back-reaction: 



Oct = / d'^xWix) -i^^en + 4A(l+(5, 



J vol 1 



-9 



(14) 



Here W{x) is their "window function" which is supposed to be a scalar that 
does not depend upon the metric, [f] Combining equations (38) and (41) of 
Garriga and Tanaka results in their definition for Rren'- 



Rr,„ = R- 16nG X ^ + 4A(5voi 



(15) 



Substituting ( ITSll in (HM reveals the key operator of Garriga and Tanaka as 
the smeared trace of the renormalized field equations: 



GT 



167rG' X / d*xW{x)g^"'{x 



6{S+AS) 
' 6g^"'{x) 



(16) 



Of course this does vanish, as do all components of the renormalized field 
equations, in all expectation values, with or without the smearing. 



• Consistency With Screening: Garriga and Tanaka have not shown their 
definition of the renormalized Ricci scalar is constant but rather that it, times 
the volume element, equals another operator: 



-9 



AA((l+6. 



vol 



-9 



(17) 



This is completely consistent with the secular screening mechanism described 
in Section 1. It is incorrect to think of R and as possessing distinct time 
dependences in flTTl) but one can do so at leading logarithm order. In that 
case, working in conformal coordinates, the mechanism of Section 1 implies 
the following relaxations for the Ricci scalar and volume element: 



i? = 4A + 87rG'(p - 3j9) 



4A 



l-4^G^H^\n{a)+0{G^) 
l-i^G^HHn\a) + 0{G^\ 



(18) 
(19) 



""We remark that W{x) must involve the metric if it is not constant. Garriga and Tanaka 
seem to feel that metric dependence in W{x) would prevent the equations of motion from 
holding, but any obstruction of this sort can be absorbed into an operator ordering [35j . 
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At leading logarithm order, the change in R is insignificant compared with 
that of the volume element and both sides of equation (ITTI) are equal: 




4A((1+Mv^), 

4Aa^{l-i^[GHHn{a)f + 0{G^)} . (20) 



This continues to be true at all orders in the loop expansion, as per relation 
(IHl), because R involves differentiated metrics whereas y/—g does not. Of 
course Garriga and Tanaka want equation flTTl) to be exact and not just valid 
at leading logarithm order. However, the necessary sub-leading logarithms 
can easily be supplied by cross terms between R and y/—g, and by the oper- 
ator 5vol- 

• Peculiar Renormalization: Although the renormalization scheme (1151) 
is consistent with secular slowing, we still find it dubious. As Garriga and 
Tanaka point out, using the equations of motion inside a functional integral 
requires the integrand to be an invariant [35]. Hence one is not considering 
"i?" and "1" but rather ^^Ry/—g " and "v^— Renormalizing such compos- 
ite operators requires that their insertion in a IPI diagram be accompanied 
by the insertion of counter-operators: 



We know of no previous study of these operators with a non-zero cosmolog- 
ical constant but, on general grounds of symmetry and the degree of 1-loop 
divergences one expects [2S1 EZ] : 



5R = IGttG X [(3,nR + f32R^ + (^^R^'R^, + [^aC"''"' C p^^, + 471 Ai?) 



The trace of the renormalized field equations ( JT2l) requires that the divergent 
part of (3i should agree with the divergent part of — 6q;i, and that (32-i should 
vanish. That is no problem. However, the renormalization scheme (1151) of 
Garriga and Tanaka also requires that the finite part of (3i should agree with 
the finite part of — 6ai. An escalating series of similar relations must hold 
as one moves up the loop expansion. 




(21) 
(22) 




(23) 



+0{G^) . 



(24) 



9 



There are two ways of viewing this. Either Garriga and Tanaka have 
solved the problem of quantum gravity by uniquely specifying the finite parts 
of every a, except a2l§ or else they are defining the Ricci scalar so as to absorb 
known classical effects. In the latter case one could also absorb the crucial 
E? term of Starobinskii's model [31] to conclude that the "renormalized Ricci 
scalar" is constant. Indeed, there seems no reason to restrict the procedure 
to pure gravity. Were it applied to scalar-driven inflation — by including the 
scalar action in the "AS"' they use in (fT5|) — then one would conclude that 
scalar-driven inflation never ends. 

It might be objected that a key distinction between the counterterms of 
quantum gravity and the inflaton potential is that the former contain factors 
of h whereas the latter does not. If so, then suppose the inflaton potential 
derives from a one loop correction, such as the Coleman- Weinberg potential 
of New Inflation [38]. The inflaton potential now carries a factor of so it 
can be absorbed into the "renormalized Ricci scalar" and we have just shown 
that New Inflation is pure de Sitter. 

• The Composite Operator R^^n Of l[T5\) Is Neither Properly Defined Nor 
Is It Constant: Correctly deflned composite operators have the property that 
inserting them in IPI functions gives flnite results. With the conventions of 
Garriga and Tanaka, the composite operators Rren\/—g and (l + 5ren)v^— ^ 
have this property but Rj-en does not. Simply divide out the volume element 
from both sides of their key relation ([9]): 

Rren = 4A(1 + 5vol) • (25) 

Now note two crucial features of the volume counterterm 6voi'- 

1. It isn't flnite; and 

2. It is an operator rather than a constant. 

The flrst feature can be seen by expanding the volume element operator in 
powers of the graviton fleld: 

V^g = a^{l + + y (n)^ - ^h^-h,. + ...). (26) 
^The single parameter a2 would not prevent us from making sense of quantum gravity. 



Rn 



4A(1 + 6, 
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Inserting ^—g in the vacuum amphtude gives a divergence at order = 
IGttG from the coincident propagator. From expression fl23p we see that this 
divergence is canceled by the 1-loop contribution to 5voi, which is k^jiR = 
'-fiD{D — 1)k,'^H'^ + O(k^). It follows that 71 is divergent, so inserting just 
(5voi into a IPI function, without the factor of ^—g, does not produce finite 
results. We therefore conclude that Garriga and Tanaka cannot claim to have 
shown that the Ricci scalar is constant — or anything else — because they 
have not correctly defined it as a composite operator. Further, the quantity 
i?ren, is uot coustaut but rather an infinite series of time-dependent operators. 

4 Ricci Scalar Constancy and Expansion 

It is time to critically examine the unstated assumption of Garriga and 
Tanaka that proving the Ricci scalar is constant implies de Sitter expansion. 
This is quite incorrect. Neither the classical field equation, R^^ — \Rgp,u + 
^9puv = 0, nor its trace, — i? + 4A = 0, precludes the presence of gravitational 
radiation. In fact, all gravitational wave solutions with A 7^ must obey 
these equations. Because gravity is a nonlinear theory, gravitational radia- 
tion interacts with itself. Because the gravitational force is attractive, this 
self-interaction opposes the expansion of spacetime. If enough gravitational 
radiation is present, the result is not expansion but rather collapse into one 
or more black holes. This is an unavoidable consequence of the singularity 
theorems. Therefore, one cannot conclude that constant Ricci scalar implies 
a constant expansion rate because this is not even true classically. 

It might be objected that a positive cosmological constant opposes the 
tendency towards collapse of too much gravitational radiation. So the cor- 
rect result is that a fixed, initial distribution of gravitational radiation may 
well collapse to one or more black holes, but the late time geometry away 
from these black holes will approach the local de Sitter form for which 
-^''(TAii/ — H^i^^9(ru — Su9(Tfj.)- That would be true enough classically, where no 
new gravitational radiation can be generated beyond what was present ini- 
tially, but it is not correct when quantum effects are included. The simplest 
way to understand quantum effects is as the classical response to the source 
provided by uncertainty principle. In these terms, our screening mechanism 
is driven by the steady injection of gravitational radiation, throughout space, 

^The divergence would be "4^ -p34 in our gauge, implying 71 = ~4gi^"D34 + finite. 
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as more and more infrared virtual gravitons are ripped out of the vacuum. 
These gravitons must slow the expansion of spacetime because that is what 
a classical distribution of gravitational radiation would do. However, such a 
classical distribution would obey R = 4A. 

5 Epilogue 

It seems obvious that gravitons contribute to vacuum energy. At 1-loop or- 
der, where each mode contributes separately, infrared gravitons participate 
on an equal footing with ultraviolet ones; this can be shown by explicit com- 
putation [HEIE]- Higher loop contributions involve non-linear combinations 
of gravitons which are best viewed as integrals over position space. These 
integrals reach back from the point of observation, along the past light-cone 
to the fixed initial value surface. Because the past light-cone grows as the 
time of observation evolves, these higher loop contributions can grow as well. 
That has to be expected in view of the infrared singularities of quantum 
gravity on a locally de Sitter background . 

Nothing ought to seem dubious about this. Although the effect has not 
yet been demonstrated for an invariant measure of the quantum gravitational 
back-reaction on inflation, the same phenomenon injects secular dependence 
into other invariant quantities p!Sl,fT9j. It also engenders secular contributions 
to the vacuum energies in scalar models [221 ESI El] for which there is no gauge 
issue to prevent one from drawing obvious conclusions. 

We have shown that the renormalization scheme employed by Garriga 
and Tanaka is highly dubious and that, even if accepted, it does not preclude 
secular slowing. Nor would proving any plausibly defined "Ricci scalar" to 
be constant preclude secular slowing. Of course the burden of proving that 
secular slowing does occur rests quite properly with us. And even an ex- 
plicit perturbative demonstration of this at 2-loop order would not establish 
quantum gravitational back-reaction either as an explanation for why the 
observed cosmological constant is so small or as the basis of a viable model 
of inflation. For that one would need a reliable way of computing in the late 
time regime, after back-reaction has become non-perturbatively strong. 

We do not think such computations are beyond reach. Starobinskii has 
proposed a simple stochastic formalism [IQ] that reproduces the leading in- 
frared logarithms of scalar potential models to all orders [2l| and which can 
be summed to give non-perturbative results [23] • We have recently extended 
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Starobinskii's formalism to massless, minimally coupled scalars that interact 
with other kinds of fields. This has led to explicit, non-perturbative results 
for Yukawa theory [29j and for scalar QED [26]. We do not yet know how 
to treat derivative interactions of the sort quantum gravity possesses but the 
problem does not seem insolvable. 
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